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(— ^ ' We study the effective interactions between a test charge Q and a one-component plasma, i.e. 

rvj I a complex made up of mobile point particles with charge q, and a uniform oppositely charged 

background. The background has the form of a flat disk, in which the mobile charges can move. 

The test particle is approached perpendicularly to the disk, along its axis of symmetry. All particles 

interact by a logarithmic potential. The long and short distance features of the effective potential - 

the free energy of the system for a given distance between Q and the disk- are worked out analytically 

in detail. They crucially depend on the sign of Q/q, and on the global charge borne by the discotic 

complex, that can vanish. While most results are obtained at the intermediate coupling F = I3q'^ = 2 

(/3 being the inverse temperature), we have also investigated situations with stronger couplings: T = 

ri ■ 4 and 6. We have found that at large distances, the sign of the effective force reflects subtle details 

O ' of the charge distribution on the disk, whereas at short distances, polarization effects invariably lead 

D to effective attractions. 
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^ , I. INTRODUCTION 



Electric charges are ubiquitous in the colloidal domain, and often a major player shaping the behaviour of soft matter 

systems. Counter-intuitive phenomena often ensue, such as overcharging (charge inversion) or effective attraction 

i-rt ■ between like-charged macro-ions [1-7]. To rationalize such observations that are the fingerprints of correlation effects, 

^ ■ simplified models are welcome, that should furthermore be treated beyond the mean- field level [8]. Interestingly, the 

O ■ physics of strongly coupled charged systems has witnessed relevant progress in the last 15 years [4, 9-13], while the 

, ^, ' study in the weak coupling limit where mean- field arguments hold, started about 100 years ago [14, 15]. The study 

[ of intermediate Coulombic couplings, though, appears more elusive [16-19] and will be the focus of our interest in the 

CS| ■ present paper. 

^ [ The system under scrutiny here is a variant of Thomson's plum pudding model (see [20-22] and references therein), 

25 ■ Eilso referred to as the One Component Plasma [4, 23, 24]. Point particles with charge q are embedded in a two- 
dimensional flat disk V of radius R. In addition, a uniformly charged background is present in the disk region (see 
Fig. 1). While the charged background is fixed, the particles are free to move in T). They interact through a log 
potential, the form taken by Coulomb law in two dimensions. The relevant coupling parameter is F = /Sq'^, where /? is 
the inverse temperature. At small F (formally F — >■ 0), the Poisson-Boltzmann mean- field description holds [25] [26]. 




FIG. 1. Definition of the system under consideration: a disk D with fixed and uniform background charge (hatched area), in 
which A'' mobile oppositely charged particles, shown by the bullets, are free to move {N = 4 in the figure). The total charge 
of the background is —Nbq, while each mobile ion bears a charge q. The total charge on the disk (background plus free ions) 
is therefore {N — Ni,)q. A test ion with charge Q approaches the disk along the symmetry axis shown, defining x-coordinate 
{x = when the test charge lies on the disk, at the center.) 

As is often the case for 2D Coulombic problems, the coupling parameter F = 2 lends itself to an exact analytical 
treatment, see Refs [24, 27, 28]. The goal is here to extend the exact analysis at F = 2 to investigate the interactions 



between the disk bearing the mobile charges, and a test charge Q that is approached perpendicularly to the disk, 
along the axis of symmetry (see Fig. 1). We shall assume that Q and all other charges (mobile + background) interact 
through a log potential. Since Q explores a third additional direction compared to those in which the point g-charges 
and the background disk are confined, the choice of such a potential can be questioned: it does not correspond to the 
solution of Poisson's equation in three dimensions. This is however the price for obtaining analytical results, that shed 
light on phenomena at work in more realistic systems. In particular, we will be interested in the effective interactions 
between the fixed charge Q and other charges, that can be seen as mimicking a colloid (the uniform background), 
dressed by a double-layer of counterions (the mobile g-charges) . 

The model will be defined in section II, where the theoretical tools will also be introduced. In section III, the 
features of the effective potential at large distances will be addressed. While most of the present analysis pertains 
to the r = 2 case, other couplings will be addressed (namely F = 4 and 6, corresponding to smaller temperatures). 
Then, the emphasis will be in section IV on short range correlations that, through polarization effects, rule the short 
distance behaviour of the effective potential. Conclusions will be drawn in section V. 



II. MODEL AND GENERAL FORMALISM 

The system is a one-component plasma [23, 24, 27, 29] on a disk V of radius R with N mobile point charges q, 
and a fixed background charge density pb — —qrn,. The system can be globally charged, since Nj, — irR^ni, can be 
different from N. The charged particles interact with the two-dimensional logarithmic Coulomb potential. 






z;,(r„r,) = -ln^V^, (1) 



for two particles located at r^ and Tj on the disk (i is an arbitrary length) . The interaction potential between a charge 
q located at r and the background consequently reads 



Vb{r) = / qpbVc{r,r')(fr' 
Jv 



TTUbq^ 



r2-i?2ff-21n| 



(2) 



where r = |r|. The disk is in thermal equilibrium with a heat bath at an inverse temperature /? = 1/kBT . We consider 
now that a particle with charge Q approaches the disk from the axis normal to the disk that passes through its center, 
see figure 1. The charge Q is held fixed at a distance x from the disk. The interaction potential between this charged 
particle and a charge from the disk located at r, will also be taken logarithmic: 

VQ{x,r) = -\n J . (3) 

It will be convenient to use rescaled lengths r = y'lrubr, x — ^irubx, i = ^Trnb£, etc. With such a choice, the rescaled 
disk radius is i? = \/Nb. The interaction potential between the background and the approaching particle is 

yQb{x)^ / QpbVQ{x,r)(fr 
Jv 

^^[{Nb+^) HNb + x^)-Nb-d^\n ^] , (4) 



where we have chosen the arbitrary constant £ such that 1 = 1. 



A. The special coupling F = 2 

Since the one-component plasma on the disk is two-dimensional with log-potential, one can use the special techniques 
developed for two-dimensional Coulomb systems [24, 27] and random matrices [30] to compute exactly the effective 
interaction potential between the disk and the approaching charge, for a special value of the Coulomb coupling 
F = /3g2 = 2. 



Let Ti be the position of the i-th particle on the disk, in polar coordinates r^ — {ri,ipi). It is convenient to define 
Zi = r^e*'''' and Zj = riC^^^ . The total potential energy of the system can be written, up to an irrelevant constant 



N 2 ^ 

l<i<j<N 



i=l 



(5) 



i=l 



where the first two terms on the right hand side account for the test charge - mobile charge and test charge - background 
interactions respectively, while the last two terms are for the mobile charge - background and mobile charge - mobile 
charge energies. When fiq^ ~ 2, up to a multiplicative constant, the Boltzmann factor reads 



N 



e-^" = e 



-/3V<5.(.)-Qg-2H„^(.,..)-.^^ -Q |^_^^ 



The product rii<i<7"<Af('^« ^ -^j) '^^ ^ Vandcrmondc determinant dct(z^ ). Defining 



1=1 l<i<j<N 



the Boltzmann factor can be written as 



g-/3ff ^ ^-PVq,{x) 



det(V'j_i(rj))i<^^^-<jv 



The functions ipj are orthogonal 



2P,ir)M^)d'r^0 ifj^fc. 



with norm 



1 
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(6) 



(7) 



(8) 



(9) 



(10) 



If Q/q is a positive integer, this can be expressed in terms of incomplete gamma functions "f{k, N^) — Jq '' t'^ ^e * dt. 
For instance, when Q = q, 



Ub 



(11) 



The configurational canonical partition function is 

TV 



Nl J-r,N JJ- 



m 



,-fSVcMx) 
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N 



dct(V'j_i(ri))^ 



<i,j<N\ 



Ud^ 



(12) 



If the determinant is explicitly expanded, and the integrals performed, the result simplifies [27, 30], due to the 
orthogonality of the functions ipj 



N-l 



(13) 



Up to an additive constant, the effective interaction potential, t4ff(x), between the disk and the approaching charge 
Q, is given by [31] e^'^^^^^'^^) oc Z, and more specifically, we choose 



-pV^ttix) 



Zo 



(14) 



where Zg is the x-independent partition function when Q ~ 0. The above definition ensures that for A^ = TVf,, 
Vcs{x) -^ when a; — >■ oo. On the other hand, for N y^ Ni,, Vi2g(x) diverges for x — >■ oo, see below. The physical 



meaning of Vcs is clear : —dVcff{x)/dx provides the mean force experienced by Q, averaged over all possible fluctuations 
of charge configurations on the disk. The function l^g- is precisely the free energy of the system, for a given test charge 
- disk distance x. Therefore, 



Q 

q 



^T4ff(x) = -^ [{Nb + x^) HNb + x^) -Nb-x^ Ini^ 



N 

E 



In / t^-\x^ + tf/''e-*dt-\n-f{j,Nb) 





(15) 



In the special case Q — q, wc obtain 



l3Vcg{x) = {Nb + x^) HNb + x^) -Nb-x^ In x^ 






X 



lU.Nb) 



(16) 



The density profile n(r) on the disk can also be obtained explicitly [24, 30], and will be discussed in some detail 
below 



N~l 



<r) = E 



j=0 II ^-J I 
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^J(j2^p2)Q/ 



^ L t^(^^+t)^^'^e~tdt 



(17) 



It can be checked that the two situations where Q — and x — > oo are equivalent, since both decouple the test charge 
from those on the disk. 



B. Arbitrary even coupling parameters 



For couplings parameters F — jSq^ — 27, with 7 an integer, the partition function of the system, and the effective 
potential, can be computed for small enough number of particles N , by using a method developed in [32, 33], based 
on techniques used in the study of the quantum Hall effect [34-36]. We provide here some details on the methods. 

Up to a multiplicative constant, the Boltzmann factor of the system reads 



^-PH ^ ^-fiVQ,(x) 



det(V'i-i(r»))^<^_^.<jY 



27 



(18) 



where, now, the orthogonal functions ipk are 



V'fc(r) = [u;(r)]i/2ife 



(19) 



with 



w(r) 



— p-^l^VQix.r)-^?^ 



(20) 



The key idea to compute the partition function is to expand [det(z^ )]'' in terms of appropriate orthogonal poly- 
nomials [33]. For 7 even, the expansion is in terms of symmetric monomials, whereas for 7 odd, it is expanded 
in terms of antisymmetric polynomials. The coefficients of the expansion are conveniently indexed by a partition 

M = (MIi ■ • • : Miv) of -fN{N — l)/2, for example for 7 even. 



[det(5r')r = E^A.Sym(zr...z^«) 



(21) 



with the symmetric monomial 



Sym(zr...zr)=n^E<(i 



) ■ • • ^aiN) 



(22) 



o-gSn 



where Sn is the permutation group of N elements and rrii is the multiphcity of the integer i in the partition /i. A 
similar expression is used for 7 odd with antisymmetrized monomials. 

Due to the orthogonality of the (anti)symmetric monomials, the partition function is finally given also as an 
expansion similar to the one of the power 7 of the Vandermonde determinant, see [33] for details. The final expression 
for the effective potential is 

pV,six)=pVQbix)-ln^ (23) 



zs 



with 



^*-En^nil^.J^ (24) 



fc=l 



1 Z"^'' 
{ry^dr = — / e-^*{x^ + ty^'H^ dt , (25) 



^A? , ,, 

Jv ^b Jo 

and Zq is Z* evaluated when Q — 0. In the case when 7 is odd, the factor Y[i "nii- — 1, since due to the antisymmetry, 
the admitted partitions /i do not have repeated numbers. 

This method can equivalently be formulated by transforming the classical problem of the one-component plasma in 
a quantum problem of a linear chain of interacting fermions, as explained in Refs. [28, 37]. The starting point for this 
method is to write the Vandermonde determinant as a Gaussian integral over Grassmann variables. The final result 
is again (23). 

For the present work, we did some calculations up to iV = 11 particles. The coefficients c^ needed for the numerical 
calculations where kindly provided by L. Samaj for 7 = 2uptoiV=10 and for 7 = 3 up to A^ = 9. For 7 = 2 
and A = 11, and 7 = 3 and A^ = 10 and A^ = 11, we obtained the coefficients using the algorithm recently proposed 
by Bernevig and Regnault [38], and their Jack polynomial generator online code [39]. We now turn to the results 
obtained from the previous analysis, starting with the effective potential experienced by the test charge Q at large 
distances from the disk. 

III. LONG DISTANCE BEHAVIOR 

A. General results at arbitrary couplings 

As mentioned earlier, the effective interaction potential, also known as potential of mean force, Vcs{x)^ has the 
property that — VT4ff is the mean force experienced by Q. It is interesting to introduce another quantity, V{x)^ the 
electric potential created by the average charge density distribution q{n{r) — nt) at the position of the charge Q 

F(x) = / q{n{r) -nb)vQ{x,r)dr (26) 

Jv 

Because of the fluctuations and the fact that the presence of Q at position x modifles the density on the disk, in 
general Vcs{x) 7^ QV{x), only for a small infinitesimal charge Q the equality holds. For arbitrary Q, a simple relation 
can be found between the two, by noticing that the total potential energy of the system (5) depends linearly on Q, 
then 

QQ =^P j q{n{r;ri,...,rN)-nb)vQ{x,r)dr e^i^" (27) 

where n(r; ri, . . . , r^r) = X]i=i ^i^ ^ ^i) '^^ ^hc microscopic density. Averaging this relation over all the configurations 
of the ions on the disk, we find 

dVMx) 



dQ 



= V{x) . (28) 



At large distances from the disk, expanding vq{x, r) for r <C a;, one can obtain the multipolar expansion of the electric 
potential V 

V{x)^-q{N-Nb)\nx-q^ + 0{l/x^) (29) 



where the relevant quadrupole moment Q2 results from the second moment of the excess density [n" (r) — nt,] 



^2 = j An\r)-nk]dv 

JV 



(30) 



where n^{r) is the density when x ^ 00 (or equivalently (5 = 0). Since we are computing the potential on the x-axis, 
no dipolar contribution remains, while the logarithmic monopole contribution stems from the global charge of the 
disk, q{N — iVh). Since up to terms of higher power than 1/x^, Eq. (29) shows that V{x) does not depend on Q, 
integrating Eq. (28) one finds that Eq. (29) also gives the large x expansion of VcS (multiplied by Q). 

We also mention here a sum rule that turns out interesting for the following discussion. The quadrupolar moment 
Q2 can be shown to be related to the mobile particle density at contact through [33, 40] 



r 



rib 



Ub 1 



n'>{R) 



(31) 



for a neutral disk. 

For the next order of the expansion, in 1/x^, the situation is more involved. The expansion of V{x) cannot be 
obtained only from the next multipole Q4 ~ L, r^{n'^{r) — Ub) dr, since the density n{r) itself depends also on x and 
one needs to take into account the next-to- leading order of the large-x expansion of n{r) to compute properly the 
expansion of V{x) up to order 1/x^. This next-to-leading order is of order 1/x^ and it is proportional to Q as it can 
be checked by expanding e^^^ for large x. It has the form 



^ N 0/ X , (3Qqd2{r) , ,^4, 

n(r) = n (r) -\ h 0[l/x ) 



(32) 



where d2{r) is a function only of r and Pq'^. Using this expansion one can obtain from Eq. (26) the expansion of V{x) 
up to the order 1/x^ 



V{x) == -q{N - Nb) Inx - q— + <l[^ ^ 



f r^d2{r)dr]^ + 0{l/x''). 
Jv J X 



Then integrating with respect to Q one finds 

V,s{x) = -qQ{N - Nb) \nx- qQ^ + qQ 



^'^^ I' r'd2{r)dr]^ + 0{l/x^). 



V 



(33) 



(34) 



The term involving the second moment of ^2 differs by a factor Q/2 between V{x) and t4ff(x). In the following section 
we illustrate these considerations on the explicit results obtained when F — j3q^ = 2. 



Results at r = 2 



Unless otherwise specified, the results reported correspond to F = 2. For Q arbitrary, the long distance behavior of 
the effective interaction (15), for N and Nb fixed, x^ ^ N and x^ ^ Nb is 



V,s{x) = -Qq{N -Nb)lnx + ^^ — 



Qql_ 

2 P 



N 



'4-Y. 



J=l 



7(j + l,^b) 
l{],Nb) 



(35) 
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E 



{Q-q) 



7(j-f2,JVb) ^ l{j + l.Nbf 
7(j, Nb) 



l{j, Nb 



+ 0(l/2?6), 



the structure of which deserves some comments. Up to order 1/x^, such a series has the form of a multipolar expansion, 
in agreement with the discussion from the previous section. Indeed, the coefficient of 1/x^ is precisely —qQQ2/2 as it 
can be checked by computing the second moment of the excess density from the explicit expression (17) when Q = 0. 
Eq. (35) can be compared to the large-x expansion of the electric potential 



V{x) = -q{N - Nb)\nx + ^ — 



2J2 



N 



^-E 



i=i 



lU + l,Nb) 
7(j, Nb) 



(36) 
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FIG. 2. Reduced charge density profile in the disk, for different neutral situations (A^ = A^;,), and x — > oo. The arrow on the 
right hand side indicates the limiting value ln2 ~ 0.693 that is reached in the large A'^ limit. The total charge density profile is 
q[n{r) — n;,]. For large A'^, it thus vanishes except in a small region of linear size l/v'iVt in the vicinity of the boundary r = R. 



where one can explicitly check that dgVcs — V. 

Let us discuss further the expansion of V^ff up to order 1/x^. Using the properties of the incomplete gamma 
function, that allow us to write the coefficient of the 1/x'^ term appearing in Eq. (35) as 



N 



^-S 



i=i 



7(j + l,JV,) _ N^-N^^ 



N ^ 



N^ 



2 Ub 



(37) 



where n°(i?) is the density of particles at the edge of the disk in the absence of the charge Q, i.e. Eq. (17) with Q—0 
aXr^R^ VTVfc. Thus, 



V^ff (x) = -Qq{N -Nb)\nx+^^ 



Qql_ 
2 x2 



N^ - N^ N Nbn^^^R) 



rib 



+ 0{l/x^). 



For a neutral disk, N ~ Nb, and 



V^cff(x) 



QqNb 1 



rib 



+ 0(1/2?*). 



(38) 



(39) 



This result is an explicit check, at F == 2, of the multipolar expansion (29) combined with the sum rule (31). One notice 
that the Vcff is repulsive for Qq > 0. This can be understood as follows. When x — >■ oo, the charge density profile 
n{r) inside the disk is the same one as for a disk alone (without the approaching charge Q), found in Refs. [24, 41]. 
The density n(r) is equal to the background density Ub in the bulk of the disk (local neutrality in the bulk). Close 
to the boundary, it raises above the background density, then falls below it [41], see figure 2. Therefore, there are 
two concentric layers of charges close to the edge: the inner layer bears a net charge which is of the same sign as q 
[i.e. n{r) > n^], while the outer one is opposite, by electroneutrality This ensures that Q2 is generically negative, so 
that the quadrupolar term yields an effective interaction (disk-test charge) that is of the same sign as Qq, i.e. repulsive 
for Qq > 0. 

When N — > 00, more explicit results can be obtained. In this limit, the density at the edge of the disk takes a 
simple form [41], n^'^^{R) = nhln2, so that 



Kff(x) 



QqNb 1 



ln2- 



1 



0{l/x^ 



QqR^ 1 



ln2- 



1 



0{l/x^ 



(40) 



with ln2 — i ~ 0.19 > 0, which is consistent with the generic discussion above on the negative sign of Q2- Again, the 
effective potential is attractive at large distances for Qq < 0, repulsive for Qq > 0, for a neutral disk. 

The quadrupolar route allows us to obtain results for strongly coupled systems (large F), making use of the sum rule 
(31). We note in passing that this general result is compatible with the value Q2 = —qR^{\og2 — 1/2) that holds at 



r = 2 when the number of mobile charges on the disk becomes large, see Eqs. (40) and (29). When F itself turns large, 
the system crystallizes, but the sum rule (31) remains valid, provided n{R) is replaced by the average of the contact 
density over the perimeter of the disk [33, 40]. It is physically reasonable to suppose that the average of n{R)/nb 
remains bounded in this limit. Actually, for the three-dimensional analogue of this model, with 1/r interaction, this 
is the case [42]. Then Eq. (31) becomes Q2 ~ —R^/2. We consequently have 

V.n{x)^^ ^\+0{l/x% (41) 

which is again repulsive for Qq > 0, attractive otherwise. It can be mentioned here that the scaling result Q2 fx —R^ 
is readily recovered by the two concentric layers simplified viewpoint. For large N , there exists a outermost corona 
void of charges: particles are depleted there, as they are in the plum pudding model, see e.g. [21, 22, 42], the width 
of which is given by the typical distance between particles 5 oc R/^/N (at F = 2, 6 is already the typical distance 
between the density maximum and the disk radius that can be seen in Fig. 2 for large N). The charge in this corona 
is given by —qribRS, which contributes a quantity —Rnj^SR^ to the quadrupole moment. This charge is compensated 
by an oppositely charged ring, located at _R — (5, which contributes a quantity R5rn,{R — 5)"^ to Q2- Summing both 
contributions, assuming that the particles and background with r < R — 5 do not contribute to Q2, and remembering 
that (5 ^ i?, we arrive at Q2 ex —ni,R?5^ ~ ~R^. A very similar argument holds at F = 2, since the two corona 
approach is also valid. 

IV. SHORT SCALE FEATURES 

We now turn to the study of the phenomenology at shorter distances, which is different depending on whether the 
particle approaching the disk has a charge of the same sign of the mobile particles on the disk [Q/q > 0), or a charge 
of opposite sign. In addition, the cases of globally neutral or charged disks should be treated separately, and the 
different cases are ruled by different sorts of polarization effects. 

A. Case Q/q > 

1. Neutral disk 

We consider a globally neutral disk N ~ N^. We study in this section if it is possible to overcharge this object, 
by approaching a particle that has a charge Q with the same sign of the mobile counterions on the disk. At large 
distances, we know from (38) that the interaction is repulsive. We anticipate that this behavior changes when the 
charge Q is close enough to the disk, since the intruder Q should then create a correlation hole, pushing mobile charges 
closer to the boundary r = R, and thereby gaining Coulombic energy from hole opened. This is the mechanism behind 
charge inversion in colloidal systems, that has been reviewed for situations of strong coupling in Ref. [3]. 

The short-distance behavior of the effective potential (15) is, when Q/q > 0, 



l3V,s{x) 
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Nb\nNb-Nb + x^ 1 + ln^ 
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Q^2 y li^^ _V +o(x\ 2?2(i+Q/,)A ^ (42) 



X' 



1 U ^[j + ^,m) 



which is clearly an increasing function of x when ir <C 1. Therefore, at close distance from the disk, the interaction 
turns out to be attractive. This can be observed in figure 3, where the effective potential indeed increases at short 
distances, reaches a maximum and then decreases upon increasing the distance between the test particle and the disk. 
Note that at x = 0, the effective potential takes a finite value 

n ^ -y(j + 0. Nh) 

/3T/eff<0) = ^ im In Nb -Nb)-Y^ In ' ,.%, , (43) 

9 p^ 7(j, Nb) 

although it is not shown in all figures. 
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FIG. 3. The effective interaction between the disk and an approaching ion with charge Q — q (left graph) or Q = IQq (right 
graph). The disk is globaUy neutral with Nb — N ions of charge q. The distance x is expressed in reduced units (x), in which 
the disk radius reads VNb, and therefore takes different values for the three curves shown. 

Let X* be the distance at which the interaction potential reaches its maximum, and x* — yJmibX* . x* is the 
minimum distance that one has to approach the charged particle in order to overcome the natural repulsion of the 
disk. The corresponding (free) energy cost to overcharge the disk is given by F^ — t4ff(x*), and V* — V'^ — T4ff(0) 
is the binding energy, i.e. the necessary energy to unbind the charged particle from the disk, once it has been 
overcharged. More generally, V* can be defined as the energy to overcome to peel off an ion from the disk. In all 
the present discussion, the energy costs alluded to correspond to the work an external operator holding the intruder 
should perform, and this equals the corresponding free energy variation of the system as a whole. 

Figure 4 shows how x* /R, V'' and V* depend on N for fixed Q/q, and on Q/q at fixed A^, respectively. First of all, 
it appears that the binding energy V* is several orders of magnitude larger than the energy cost V' . This means that 
V* ~ |Vi3ff(0)|. Second, the threshold distance x* scales like R, when N becomes large enough, a fact that is masked 
in Fig. 3 by the choice of units made (tilde variables, for which R = \fNl). More precisely, from the numerical data 
of Fig. 4, we explored how x* / R depends on the charge Q. We found, numerically, the approximate relation 



x 




+ 0(1/V^) 



(44) 



with a = 1.5 + 0{l/\/N). The fact that a is of order one means that the effect of effective attractions holds up 
to rather large distances, on the order of the disk radius. Another feature visible on Fig. 4 is that the energy cost 
increases with Q/g, but quickly saturates to a finite value V^'Jf*- On the other hand the binding energy increases as 
the charge Q increases, as expected, but it also increases with the number of mobile ions on the disk N . For N — > oo, 
using Stirling formula for the incomplete gamma functions in Eq. (43), one can obtain the analytical behavior of 
Vcs{0): and therefore the one of the binding energy, remembering that V* ~ |V^ff(0)|. We find 



/3Feff(0) 



\tlN + 0{l). 



(45) 



Figure 5 shows how the previous quantities behave under different couplings (T — Pq^ — 2,4,6). The qualitative 
features appear to be robust: The behavior is similar to the one when F = 2, with changes in the numerical values 
of X* , y^ and V* . As F increases, V^ and V* increase, and x* decreases slightly: the test ion has to come closer to 
the disk, and requires more energy to overcome the long distance repulsion. Once it overcharges the disk it is more 
energetically bounded to it. 



2. Charged disk 



If iV < Nb, the disk has a charge q{N — Nb) of opposite sign to that of the approaching ion Q. In this case the 
effective potential is attractive at all distances. Let us consider the more interesting case where the disk is already 
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FIG. 4. Left: The binding energy V* ^ the energy cost V^ and the distance x* to overcharge the globally neutral disk with an 
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a function of intruder charge Q, for N — Nb = 22. 




FIG. 5. Left: The binding energy V* , the energy cost V^ and the distance x* to overcharge the globally neutral disk with an 
additional particle of charge Q = lOg as a function of the number of particles N = Ni, on the disk, for different values of the 
Coulombic coupling F = /3g . Right : same as a function of Q, for a globally neutral disk with N = Nb — S. 



overcharged, with a net charge of the same sign as Q, i.e. N > Ni,. The question is to study the effective potential 
profile, and the distance range where it corresponds to an effective attraction. 

At large distances, the effective interaction between the ion and the disk is repulsive, and diverges as —Qq{N — 
Nb) \nx. However, from Eq. (42), we find that, at short distances, the effective potential becomes attractive. Therefore, 
as in the previous situation, there exists a distance x* below which the test charge will be attracted, which results 
in a further charge inversion of the disk. Figure 6 shows the effective potential in this situation for different charges 
Q and charges of the disk q{N — Nb). Here, one can also define a binding energy V* = T4ff(a;*) — ycff(O), necessary 
to pull out the ion from the disk once it has been "adsorbed" . Figure 7 shows how x* and V* depend on the charge 
Q of the ion and on the charge of the disk, respectively. It is also useful to emphasize that when the disk complex 
is not neutral, one cannot define the energy barrier y^. Indeed, this quantity was defined in the neutral case as 
the barrier to overcome to approach the test charge from a; = oo, down to the distance where attraction sets in. 
When N y^ Nb, the large distance effective potential diverges as {N — Nb) log x, which precludes the definition of 
V^ = Ves{x*) — T4ff(oo). This feature is absent in three dimensions, where charges interact though a 1/r potential 
(hence the possible definition of V'' also for non neutral complexes). As a consequence, the study of overcharging is 
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somewhat less rich in the present case than for three dimensional systems and overcharging is, with a log potential, 
necessarily a phenomenon of small amplitude (if not infinitesimal) : another way to rephrase previous remarks is that 
V^ diverges as soon as N ^ N},. 

To understand the mechanism behind the attraction at short distances, it is instructive to study the density 
distribution of particles in V, as Q approaches the disk. Figure 8 shows the density profile for different distances x. 
It can be seen that the correlation hole alluded to earlier is increasingly marked, when x becomes smaller: the mobile 
charges q of the disk feel the repulsion due to the charge Q, and move towards the edge of the disk. This results 
in a local negative charge density in the center of the disk, which is finally responsible for the attractive interaction 
between the disk and the intruder charge Q. 
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B. Case Q/q < 

We now turn to the case where the test particle and mobile ions on the disk have charges of opposite signs. 



1. Neutral disk 

For a globally neutral disk {N — Ni,), we know from Eq. (39) and the analysis of section III that the effective potential 
is attractive at large-distances. This behavior remains at short-distances, as illustrated in figure 9. Therefore, the 
neutral disk has a natural tendency to overcharge. 
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2. Charged disk 



If A^ > Nb, the disk has a charge q{N — Nb) of opposite sign as the approaching ion Q. In this case the effective 
potential is always attractive, a situation that is not of particular interest. We concentrate instead on the case where 
the disk has a net charge of the same sign as Q, i.e. Ni, > N. Due to this excess charge, the effective potential with 
a charge Q of the same sign as the disk is expected to be repulsive at large distances, see equation (38). However, as 
we shall see below, there is here also a change in the behavior of the effective interaction at short distances, where 
the force between the disk and the particle becomes attractive. 

The situation seems at first sight similar to the case studied in section IV A. There are some notable differences 
though. If Q/q < —1, the effective potential diverges when a; — > 0. Indeed, for Q/q < —1, equation (42) is no longer 
valid. The dominant contribution is given by the term j = 1 in the sum (15). Explicitly, it yields, for Q/q < —1, 



Kff (a;) -^ -qiQ + q) In X, 







(46) 



and 



Vas{x) 



— :- In In — 



0, 



(47) 



if Q = —q. In both cases, the small x behaviour is attractive. We note that the divergence of K:ff(a;) for x -> stems 
from the fact that the Boltzmann weight exp(/?Q(7logr) — l/r"'"*?/' is non integrable at r = whenever Q/q < — 2/r. 
As announced earlier, the effective potential is repulsive at large distances and attractive at short distances. We can 
thus again define the distance x* at which the potential becomes attractive if we approach the ion below x* . However, 
in the present case, the binding energy V* is infinite because Vnnx^o V'eff (a;) = — oo. Likewise, we cannot define 
the energy barrier V'' since either attraction applies at all distances (neutral complex on the disk), or the effective 
potential diverges at infinity (charged case). 
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The mechanism behind the attraction between these like-charged objects at short distance is now due to an accu- 
mulation of mobile charges near the center of the disk. This can be seen in Fig. 10, which shows the density profile 
of the mobile charges q, ie. the counterions. As the charge Q approaches, it strongly attracts the counterions to the 
disk center, and this results in an effective attraction, the precise form of which is non trivial. As seen on Fig. 10, for 
small X, there is a large density of counterions close to the center of the disk, but which is concentrated over a disk 
of radius of order 1 in r units. The precise analytical behavior of the charge density can be extracted from Eq. (17), 
for ir <C 1, taking into account that for Q/q < — 1 the first term of the sum is the dominant one 



n{r) ^ —rib e "" { hi 



1 



-Sj^ 



(48) 
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The total integrated charge of this counterion cloud close to the center plus the background turns out to be equal to 
—Q, as one might expect. However, the behaviour of the effective potential encoded in Eqs. (46) and (47) exhibits 
a different attraction than the bare Coulombic form — Q^logx, that would be obtained assuming the attracted 
counterions are located as a point charge at r = 0. Since they are spread over distances larger than x, the behavior 
of the potential, although logarithmic, turns out to have a different prefactor, see Eq. (46). 

The case when — 1 < Q/q < is somewhat different. The effective interaction potential is no longer logarithmic 
and has a finite value at a; = 0. For a; — > 0, 

~2(l + Q/q) 

PV,s{x) = /3Kff(0) + y- r + 0(x2) (49) 

[^ + l)l{l + Qlq,Nk) 

with T4ff(0) given by Eq. (43). Notice that since < 1 + Q/q < 1, the potential is again attractive at short distances. 
Also the power law a^2(i+Q/ij) jg different from the one of the case Q/q > where it was jp. In this case, —1 < Q/q < 0, 
it is again possible to define the binding energy, that diverges when the limit Q/q — >■ —1^ is approached. One can 
indeed show that /3Vcs{0) -- In (1 + Q/q). 

V. CONCLUSION AND DISCUSSION 

We have introduced a classical system that exhibits some of the phenomenology at work in more complex colloidal 
suspensions. An ensemble of N point particles with charge q are free to move within a disk of radius R, that 
bears a uniform background charge of surface density ~qNh/{'KR^). The corresponding complex (mobile charges and 
background) forms a one component plasma, with a global charge {N — Ni/)q. A test point charge Q is then approached 
to the complex, perpendicularly to the disk plane, along its axis of symmetry (x-axis, see Fig. 1). All charges were 
assumed to interact through a log potential, a choice that is convenient for the derivation of analytical results and 
for the discussion of the physical mechanisms, but that we emphasized as somewhat irrealistic for a real Coulombic 
problem in three dimensions. We have studied in detail the x-dependent effective potential VcS experienced by the 
intruder Q, defined as the free energy of the complete charge distribution for a given distance x between the test 
charge and the complex. 

At short distances x, 14ff is always attractive, with different underlying mechanisms depending on the sign of Q / q. 
If the intruder and the mobile charges are like-charged, the intruder creates its own correlation hole as it approaches 
the disk. The resulting short-range attraction resulting from this polarization is analogous to its three dimensional 
counterpart explaining charge inversion (overcharging, see Ref. [3]). If on the other hand Q/q < 0, the test particle 
attracts an excess of mobile charges in the vicinity of the disk center, which overcomes the background - test charge 
repulsion. In this case, we found a diverging attraction for Q/q < —1, which precludes the definition of a binding 
energy (cost to drag the test charge away from the disk, starting from a; = 0, the point of contact). 

The long distance behaviour is also of interest. If the complex has a net charge, the leading contribution to V^ff 
reads —Qq{N — iVb)lnx, which leads to the expected like-charge repulsion at large x. The neutral case N = Nb is 
more subtle, and it has been shown that a key quantity to rationalize Vcs is the quadrupolar moment Q2 of the total 
charge distribution on the disk. At large x, polarization effects disappear, and the mobile charges adopt a profile that 
compensates for the background charge in the bulk of the disk, while they are expelled from the immediate vicinity 
of the disk edge r — R, thereby creating a charge imbalance far from the disk center only. This necessarily leads 
to a negative value of Q2, and hence to a repulsive behaviour at large x, when Q/q > 0. Indeed, what matters for 
large distance interactions is the charges that are closest to the intruder, and they happen to be the mobile charges 
expelled from r — R (see the region where n > nt, in Fig. 2, for A^ = 50 or A^ = 500, or equivalently, see the arrow 
in Fig. 11 below). The ensuing interaction is repulsive when Q and q are of the same sign. This leads us to a final 
remark that illustrates the subtlety of the long distance effective potential. Consider a variant of the previous model, 
where the mobile charges are no longer confined in the disk r < R, but can explore the full disk plane (they are thus 
still 2D confined, but unbounded in the plane). The uniform background, as before, is a disk of radius R. We can 
repeat the analysis for F = 2, which leads to a profile n{r) that departs from the one reported above in an essential 
way: As can be seen in Fig. 11, it is monotonously decreasing, as happens to be the case at mean-field level [43] (i.e. 
for F ^ 0). For A^ > Ni,, the decay of the density profile n{r) at large distances is algebraic in 1/r^ [44, 45], leading 
to a divergent quadrupole. Furthermore, the density profile of this "unbounded" model shows a peculiarity, when 
N > Ni,, Jjj2 n{r) (Pr — Ni, — 1. Since there were originally N mobile particles, this means that A^ — A";, + 1 particles 
have escaped to infinity. This can be checked explicitly at F = 2 [44-46], but more generally, it is a manifestation of 
the Onsager-Manning-Oosawa condensation phenomenon [47-49]: only a fraction (A"f, — 1)/Ni, of the mobile ions are 
"condensed" inside or in the vicinity of the disk. This is a consequence of the logarithmic interaction between the ions 
and the disk when they are outside the disk. As a consequence, and at variance with the bounded model where the 
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charges stay in the disk, the global charge of the complex (disk + mobile ions) is —q, whenever N > Ni,. Therefore 
one expect that the effective interaction of this complex with the charge Q at large x will be attractive for Q/q > 0. 
This situation is opposite to the one met with the bounded model. 
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FIG. 11. Density profile of counterions in the disk, for the "unbounded" and "bounded" models (counter-ions are either allowed 
to explore the region r > R, or not). Here A^ = A'^t = 40 and F = 2. The counter-ion excess for the bounded model -shown 
by the arrow, and studied extensively in this paper- leads to a negative quadrupole moment Q2, see section III, while in the 
unbounded case, one mobile ion escapes to infinity, leaving the complex (disk -I- ions) with a net charge —q. This results in 
large distance effective forces on the test charge that have opposite signs. 
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